Winding expansion techniques for lattice QCD with chemical potential 
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We analytically derive a decomposition of the lattice fermion determinant for Wilson's Dirac 
operator with chemical potential into winding sectors, i.e., factors with a fixed number of quarks. 
Dividing the lattice into four domains, the determinant is factorized into terms which can be classified 
with respect to the winding number of the closed loops they consist of. The individual factors 
are expressed in terms of subdeterminants and propagators on the domains of the lattice. We 
numerically analyze properties of the factorization formula and discuss two applications for the 
determination of canonical partition functions with a fixed quark number: A speedup for the Fourier 
transformation technique through a dimensional reduction, and a power series expansion. 

PACS numbers: ll.15.Ha 
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II. PARTITION OF THE LATTICE INTO 
FOUR DOMAINS 
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The lattice formulation provides a fully gauge in- 
variant non-perturbative approach to QCD. An im- 
portant aspect of the lattice discretization is the pos- 
sibility of using numerical Monte Carlo techniques 
to obtain non-perturbative results. An application 
where the success so far was rather limited is lattice 
QCD with chemical potential. The reason is the sign 
problem of the fermion determinant with a chemi- 
cal potential, which forbids a direct application of 
the Monte Carlo method such that one is reduced to 
reweighting methods. 

In these attempts the use of canonical partition 
functions with a fixed quark (or baryon) number 
has started to play an important role recently [l|- 
[J. For their use it is necessary to have an efficient 
method for evaluating the canonical partition func- 
tion at a fixed quark number. This problem is in 
turn related to organizing the fermion determinant, 
which may be viewed as a collection of closed loops, 
with respect to the winding number of these loops 
around compactified time. 

In this article we present an explicit decomposi- 
tion of the fermion determinant for Wilson's Dirac 
operator into sectors with definite winding numbers. 
The approach is based on a partition of the lattice 
into four domains, and the contributions to the in- 
dividual winding sectors may be written in terms of 
propagators on these lattice domains. The terms in 
in our factorization formula for the fermion deter- 
minant are analyzed numerically. We discuss two 
possible applications of our results to the evaluation 
of the canonical partition function at fixed quark 
number. 



The starting point of our approach is a partition 
of a X /3 lattice A into four domains with different 
ranges of the time component x^: 

A^i) = {{x,xi)\xe As, -(3/2 + 1, ... - I}, 

A^^'> = {{x,Xi)\xeAs,Xi^O}, 

A(3) = {{x,Xi)\xeAs,Xi^l, .../3/2-1}, 

A(4) = {(f,x4)|f e A„x4-/3/2}. (1) 

We stress that /3 is a positive integer, the inverse 
temperature in lattice units. A^ denotes the spatial 
part of the lattice. The full lattice A is the union of 
the four sub-lattices A^^'^ i — 1,2, 3, 4. The partition 
of the lattice is illustrated in Fig. [TJ 

Wilson's lattice Dirac operator D{x, y) which we 
use here can be decomposed into pieces Z?^*-* acting 
in an individual domain A*-*' and terms D^^'^^ that 
connect the neighboring domains A(*',A*^^\i =^ j. 
For that purpose we need the operators 

1,2,3,4, (2) 

- -1,1/4 = 0, 

0,y4 = -1, 

0, y4 = 1, 

1, y4 = 0, 
:/3/2-l,y4 = /3/2, 
:/3/2,y4 = /3/2-l, 
:/3/2,t/4 = -/3/2+1, 
:-/3/2+l,y4 = /3/2, 

where we also display which values the arguments x 
and y may assume. The term D^^^ is given by (we 
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FIG. 1: Partition of the lattice into four domains 
A'*',i — 1,2,3,4: The vertical direction in the plot is 
time, while the dots in horizontal direction represent a 
whole time slice. On the left-hand side we indicate the 
values of the time argument. The lattice sites in the 
same domain are enclosed in a dashed contour, and we 
label the domains on the right-hand side. 



set the lattice constant to 1 throughout this paper 
and all numerical results are given in lattice units) 
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n4 = -f +2 

where we have defined the hopping parameter n = 
1/(4 + m) with m being the bare quark mass [§]. 
We use the convention U-j{x,X4) = Uj{x — j,X4)^, 
where j denotes the unit vector in j-direction. The 
term I?'^' is obtained by shifting the range of the 
auxiliary sums running over 714 in the last two terms 
of dsn by (3/2. 

The terms D'-^-' and D'-^^ live on only a single time 
slice, with Z3*^^) given by 



(4) 
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Uj{x, X4)Sg^-^ffda,^.oSy^^o. 



Z)('*^ is obtained from D'^^\ by setting the time ar- 
guments in the temporal Kronecker deltas to (3/2. 
The terms that connect the domains A^*^ have 



only hops in time direction, with 

D^'^-^\x,y) = -K^— ^4(^,-1) fej? (5:^^,-1(5^4,0, 

D^^''^'>{x,y) = - K^-^^^U4{x,-iy Ss^yS^^fiSy^^-.! , 

(5) 

and the terms D^^^^) ^ £)(3,2) ^ ^,(3,4) ^^^^ ^(4,3) g^j.g 
tained from D^^'^^ and D^^'^^ by setting the time 
arguments according to the values listed in 

The terms D^^'''' and D<^-^'^\ which periodically 
close the lattice in time direction, have a slightly 
different form, since here the chemical potential fi is 
coupled and the anti-periodic boundary conditions 
in time are implemented. We write the two terms as 

= ef^^D^^''^ , = e'^^f'Di''^^ , (6) 

with the yit-independent pieces 



Dlt^^\x,y) = K^-Y^U4{x, (3/2) Ss,ffS^^ ,3 Sy^ _^^^, 
D^o^'*\x,y) = i^^^^U4{x,(3/2)^Ss,ffS^^^^fi^^5y^,^. 



(7) 



Combining all terms D^*^ on the individual domains 
and the connecting pieces one obtains the 

standard Wilson Dirac operator with the chemical 
potential introduced as a boundary condition. 

For later use we remark that the individual terms 
obey the following 75-hermiticity properties: 

^Wt = j,D^'^j5 , i = 1,2,3,4, 

^(^.■)t ^ ^ (*,j) = (l,2),(2,3),(3,4), 



(4,1) t 



75 ^^^^'^^5 
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We remark at this point, that our approach may 
be implemented in the same way for clover improved 
Wilson fermions and also staggered fermions. For 
the latter case a dimensional reduction formula ob- 
tained in a different way can be found in fl|, Q ■ 



III. FACTORIZATION OF THE 
DETERMINANT 



We now use the partition of the lattice into do- 
mains to factorize the fermion determinant. The 
basic step is to write the Grassmann integral for the 
determinant in terms of the Grassmann fields orga- 
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nized according to the domains A^'^: 



» , 4 

'J ^ — 1 A — 1 



-V^^'^i^(3,4)^(4)_ ^(4)^(4,3)^(3) 
-^*'^I?(44)^,(1)„ ^(1)^(1,4)^(4) 



(9) 



Here we use vector /matrix notation for all indices, 
color, Dirac and space-time. The superscripts at- 
tached to the Grassmann variables ■0^*'', V*'*'' denote 
which domains of the lattice these spinors live on. 
The sum of all terms in the exponential function 
reproduces Wilson's fermion action with chemical 
potential. 

We now first integrate out the Grassmann vari- 
ables in the domains bS^^ and A(^), then in domain 
A(2) and finally in bS^^ . Solving this chain of Gaus- 
sian integrals, one may write the determinant as 



det[i:»] = AqW , 



(10) 



with 



Ao = det [i^(i)] det [I?(3)] det [Z?(2)] det {D^% 



W 



det 
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_^(4) [eA</35(4,2)_^5(4,2) 

x5(2)[e-^^5f'^)+5f*) 
We have introduced the following abbreviations: 

5(1) = (i?(l))-l , ^(3) = (i^(3))^l , 

5(4) = i?(4)-i?(4^1)5(l)i^(1^4)_ ^(4,3)^(3)^(3,4) 
g(2) ^ (5(2)) -1 ^ ^(4) ^ (5(4)) -1 ^ 
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^(4.1)^(1)^(1,2) ^ 

^(4,2) ^ ^(4,3)^(3)^(3,2) ^ 

5(2.4) ^ ^(2.1)^(1)^(1.4) ^ 

5(2.4) ^ ^(2.3)^(3)^(3,4) _ 



(12) 



The individual terms have an interesting hierarchi- 
cal structure: L)(i) and -D(^) simply are the Dirac 
operators on the lattice domains A(i) and A(^). The 
corresponding propagators are denoted by S'(i) and 
5(3). 

To leading order the operator £)(2) (and similarly 
I)(4)) is also the Dirac operator on the corresponding 
domain A(2) (and A('') respectively). However, also 
two additional correction terms appear: In the first 
term the temporal links in D^"^'^^ allow a quark to 
hop into the neighboring domain A(i). There it may 



propagate with the propagator S^^\ and in the end 
is again transported back into A(2) by the link terms 
in Z)(i'2). The second correction term is built in the 
same way, now allowing the quark to propagate in 
A(3). The term D'^^^ has exactly the same structure: 
It adds to the Dirac operator in A(*) additional terms 
that allow a quark to visit the neighboring domains 
A(i) and A(3). 

A useful way of thinking about the individual con- 
tributions is in terms of paths: The propagators 
5(1) and 5(3) may be viewed as collections of paths 
that connect any point in A(i) (A(3)) with any other 
point. Along these paths ordered products of the 
gauge links U^(x) appear, such that the individual 
contributions are gauge covariant. The correction 
terms in ^(2) and £>(") thus are paths in the neigh- 
boring domains A(i) and A(3) which are connected 
to the time slices A(2) and A("') through temporal 
hops. 

It is important to note, that the four determinants 
in the factor Aq (fTTj) do not depend on the chemi- 
cal potential /i. The only term where the chemi- 
cal potential could appear in is the first correc- 
tion term in D^'^y However, since there both Z?(4'i) 
and _D(1'4) appear, the factors exp(±/i/?) cancel, and 



only the /^-independent terms D^^'^'' and Dq''''^' re- 
main. As a consequence one expects that the four 
determinants are real. For det[D^^'>] and det[D'-^^ 
this follows from the 75-hermiticity of £>(!) and Z?(3) 
(see Eq. (|5])). Also the corrected Dirac operators 
_D(2) and Z)(4) are 75-hermitian, as can be estab- 
lished with the help of the relations in ([8]). The 
corresponding determinants thus also are real. 

We conclude, that the complex phase must be con- 
tained in the factor W defined in (fTl]) , where we ex- 
plicitly see terms that do depend on the chemical 
potential /i. In the factor W new types of matri- 
ces appear, such as the matrix Z?^'*'^-': This matrix 
connects the domains A('*) and A(2) by first hopping 
from A('i) into A(i) with the temporal hops in dI^'^\ 
subsequently transporting a quark with S^^\ and 
finally connecting to A(2) with the hops in D^^'^K 
Thus _dJ^'^^ may be considered as a large forward 
hopping term from A('*) into 

A(2) going through A(i) . 
In the same way D^'^"^ can be viewed as a large back- 
ward hopping term from 

A(4) 

into 

A(2) going through 
A(3). The other two new terms 5^ and I)f ''^^ 
have a similar interpretation as large forward and 
backward hopping terms. 

We stress at this point, that the paths which the 
large hops consist of have a length of /3/2 or longer, 
i.e., they consist of at least (3/2 individual steps. Let 
us discuss this property for the example of D^^''^"^: 
According to (fT^ this term is given by the the prod- 
uct D^'''^^5(i)D(i'2). The first and the last factor 
each contribute one temporal step. The propagator 



l(1.4) 
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5(1) is a collection of arbitrary paths, but here we 
only need to consider the paths that extend from 
the timeslice with X4 = —f3/2 + 1 to the timeslice 
at X4 = —1, since otherwise the terms D^'^^ and 
D^i'^) could not attach to S^^^ (compare Fig. [TJ. 
The paths that connect the two boundaries of A^^^ 
have at least a length of /3/2 — 2 steps, which es- 
tablishes, that o'f'^^ indeed consists of paths which 
have at least the length (3/2, and the same holds for 
T^f and i^f^). 
Finally we remark that also the large hop- 
ping terms are related to each other through 75- 
hermiticity: 



r,(2,4)t n(4.2) 



D. 



(2,4) t 



75 D'^ 75 • 
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Effectively the factor W has the structure of a 
determinant with only two timeslices given by the 
domains A(^) and A('*). With the large hops D^i'^\ 
D^^''^\ D^j^'^^ and D^'^^ we can hop between these 
timeslices. The terms 5('*) and S'(^) that appear in 
W between the large hops allow for a local propaga- 
tion within A(2) and A('*). 



IV. DECOMPOSITION INTO WINDING 
SECTORS 

The factor W of may be expanded with the 
trace- log formula for determinants. For that purpose 
we introduce the abbreviations 

Ho = 5(4)i?i'''^5(2)i?f'4)+^(4)i^(4,2)^(2)^(2,4)^ 

- S^^^Di'''^S^^^D[''^\ (14) 

These terms may all be viewed as paths that connect 
two points in A(''). In Ho these paths visit A(2) by 
either going through A(i) or A('^) and then coming 
back through the same domain. The paths in H±i 
visit all four domains, i.e, they wind once around the 
lattice. In iJ+i they wind in forward direction, while 
in H-i they run backwards. According to the above 
discussion, in all three terms Ho,H±i the paths have 
a minimum length of (3 steps. 

With the help of the definitions in and the 
trace-log formula for determinants we can now write 
W as 

W = det [l-Ho- e^'^iJ+i - e^^'^iJ-i] (15) 
= ^''^P (-E [Ho + e^('H+, + e'^'^i/^i] 

= exp (-^ e^^^ E ^ E [H,,H,,...hA. 

\ qeZ n=l ki + ...+kn=q J 



In the final step we have rewritten the n-th powers 
by introducing two auxiliary sums: The first sum 
runs over an integer variable q which below we will 
identify as the winding number of the individual loop 
contributions. The second, restricted sum runs over 
n variables ki G {—1,0, +1}. 

Due to the trace the individual terms 
Ti[Hk^Hk^...HkJ\ correspond to closed loops 
built from the paths in the i/fc. . Since if+i and 
H^i contribute a path that winds once in forward 
(backward) direction and Ho does not wind at all, 
the sum q = ki + k2 -\- ...kn is the winding number 
of the loop. It is obvious from the last expression 
in (fT5|) that the chemical potential couples to the 
winding number as expected. 

From ITSl) follows the factorization of W: 



(16) 



where the factors for the individual winding sectors 
are given by 

M/(«)= exp|-e'^«''f;i ^Tr [Ffe,fffe,...Ffc„]|. 

\ n=l fci+... + fe„=g / 

(17) 

We remark that there is a relation between the 
sectors of opposite winding number q, 



Ty(-9)(^) = iy(«) (-;/)* 



(18) 



where the asterisk denotes complex conjugation. 
This follows again from the 75-hermiticity relations 
([8]) which can be used to establish 

[Hk,Hk,...Hk„] = [Hk'Hy,^_...Hk'S ^ (19) 

where k'l — — fc„, k'2 — — fcn-i, ... k[^ — —ki. This 
result immediately leads to (fT5|) . 

For zero chemical potential the relation (fT5|) re- 
duces to 



i^(-9)(^ = o) = w^(«)(m = o)* 



(20) 



which ensures the well known fact that the determi- 
nant is real for vanishing chemical potential (VF(°) 
is real anyway). It is quite interesting to note, that 
due to the relation ((20)) at /it = the contributions 
from the non-trivial sectors with g 7^ are positive, 
since VF(-9)(0)W^(«)(0) = IW^(«)(0)|2. Thus a pos- 
sible negative sign of the /i = determinant must 
come from the non-winding contributions AoVF(°). 
For ^ the relation ([^01) is violated due to the 
different weight factors exp(±/i(7/3), and the deter- 
minant becomes complex. 

Let us finally discuss an important aspect of our 
result: We have noted above that the paths in the 
terms Hk^ consist of at least [3 steps. Thus a con- 
tribution Tr [_fffei-fffc2----fffc„] consists of closed loops 
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\W\ n 


1 Ut!L[-L-'J y^— 0.1 


\W\ n , 

1 yy Im— 0.1 




1 yy 1^—0.2 


0.05 


0.724x10'*'^ 
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3.202 
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2.743 


0.10 


0.179x10*^ 


0.509 xlO*'5 


2.838 


0.492 xlO*'^ 


2.740 


0.440x10**^ 


2.453 


0.20 


0.377x10'*'' 


0.875x10*'' 


2.319 


0.848x10*" 


2.249 


0.769x10*" 


2.039 



TABLE I: Values of Ao, |det[D]| and \W\ for three different quark masses, m = 0.05,0.1,0.2, and three different 
values of the chemical potential. At = 0, 0.1, 0.2 for a single configuration. 



which have at least a length of n (3. Long loops are 
suppressed exponentially in n, due to the fluctua- 
tions of the gauge links and the factors of n. Ob- 
viously the sum over n corresponds to a generalized 
hopping expansion. 



V. NUMERICAL ANALYSIS OF THE 
FACTORIZATION FORMULA 

Various aspects of our determinant factorization 
can be analyzed numerically. In particular we study 
the behavior of the factors and as a function 
of the mass parameter m and the chemical potential 
fi. For our exploratory analysis we use quenched 
[l0| gauge configurations on lattices of size %^ x 4. 
The configurations were generated with the Liischer- 
Weisz gauge action @ and at the value oi (3 = 7.60 
which we use here, the lattice spacing was deter- 
mined to be a = 0.193 fm. This gives rise to a 
cutoff of roughly 1 GeV and corresponds to a tem- 
perature of 254 MeV. 

As a first test we programmed both sides of the 
exact relation (fTO|) in order to test the correct imple- 
mentation of the terms in pT|) . We found agreement 
to machine precision. 

In the actual analysis of the properties of the fac- 
torization formula we begin with evaluating the fac- 
tors Aq and W (see Eq. (Illl) ) as a function of mass 
m and chemical potential /x. We refer to Aq as the 
hulk factor, while W is the winding term of the de- 
terminant. For the numerical study we consider the 
mass values m — 0.05,0.1,0.2, and scan the chemi- 
cal potential in the interval fi S [0.0, 1.0] in steps of 
0.01. Concerning the chemical potential we remark, 
that once the matrices Aq, Hq, H±i are evaluated, it 
is very cheap to compute the /z-dependent winding 
term W for many values of /z. The reason is obvi- 
ous from the definition of W in the first line of (|15p . 
which shows that the determinant one has to eval- 
uate is only over the indices of a single time slice. 
This dimensional reduction speeds up the numerical 
evaluation of the determinant considerably. 

To give an idea of the relative size of the different 
quantities, we show in Table 1 the values of bulk 
factor ^0 together with the absolute values | det[Z3]| 
and \W\ for three different masses m = 0.05, 0.1 and 
0.2, for a single gauge configuration. Note that Aq is 



independent of /i, while | det[I?]| and |iy| do depend 
on /i, and we thus give the latter two factors for 
three different values of the chemical potential, fi = 
0.0,0.1,0.2. 

What is immediately obvious from the table, is the 
rather different scale of the two factors Aq and W, 
which according to PH)) constitute the determinant, 
det[D] = AqW. While the bulk factor Aq is large 
of order 0(10^^°) - 0(10^''), depending on the mass 
parameter m, the winding term W is of order 1 (at 
least for the values of ^ in the table). Furthermore 
Aq strongly depends on the mass, while W shows 
only a weak dependence. We thus conclude that 
the real valued bulk factor Aq contributes the mass- 
dependent scale of the determinant, while W brings 
in the /i-dependence and the complex phase. Thus, 
besides separating the /i-dependence, our decompo- 
sition provides a natural factorization into the rele- 
vant scales for the bulk behavior and the response 
to the chemical potential. 

In Fig. [5] we now analyze how the winding term 
W changes with the chemical potential fx. In the top 
plot we show the absolute value \ W\ as a function of 
/i, while in the bottom plot the phase is displayed. 
Note that since Aq is real, the phase of is also 
the phase of det[_D] (up to a possible minus sign that 
could appear for exceptional configurations). The 
data in Fig. [2] are for a single configuration which 
we analyzed for three different values of the quark 
mass, m = 0.05,0.1,0.2. 

The plots show a non-trivial dependence of W on 
fi. Starting out with a weak /^-dependence for small 
fi, W changes considerably as /x gets close to the 
cutoff at /i = 1. With increasing fi the phase (j) runs 
twice through all of the interval [— tt, tt] (note that 
when exceeding tt we project the phase back into 
[—TT, tt]). What is surprising is that the mass depen- 
dence for both, \W\ and the phase cj) is rather small. 
This supports the above sketched picture of Aq car- 
rying most of the overall scale and mass dependence. 

We remark that we repeated the analysis of Fig. [2] 
for several configurations and found that the scale 
aspects of the behavior we describe in the last para- 
graph are universal: Aq carries most of the overall 
scale of the determinant and also of the mass de- 
pendence. The /.t-dependence in W shows a much 
smaller variation. Other aspects are specific for indi- 
vidual configurations, in particular how quickly and 
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FIG. 2: Absolute value (top plot) and phase (bottom) of 
the winding term W as a, function of /i for three different 
masses on a single configuration. 

with which orientation the phase ((> is run through 
as a function of fi. Also the details of the behav- 
ior of 1 14^ I as a function of ^ may vary for different 
configurations. 

Let us now analyze the buildup of the contribu- 
tions in the individual winding sectors. For that 
purpose we write the factor for winding num- 

ber q as 

1^(9) = exp (-e'^«''T(«)) , (21) 
where we define 

oo 

Til) ^ Y.Ti"^ , (22) 

ji=i 

t!,'^ = - E Tr[i7,,i/fe,...i/feJ. 
n ^ — ' 

ki + ...+kn=q 

It is straightforward to evaluate the coefficients ri'"* 
and for illustration purposes we list a few terms: 

T^^^^TtHo, (23) 
Tf ^ = TtH+iHoH_i + TtH+iH_iHo + ^Tri?-?. 



FIG. 3: The absolute values {Ti"' \ for various values of q 
as a function of n. We use two different masses, m = 0.05 
in the plot at the top, and m = 0.20 at the bottom. We 
superimpose the data from three different gauge config- 
urations. 

T^^^ = Ti-H+i , T^'^ = TtH+iHo, (24) 

' = , T^^^ ^ -Trij2^ , rj^^ = TriJ^^ffo- 

^ (25) 

Tp ^ = , T^^'' = , T^^'' = ^TriJi^i. (26) 

From q — ki + k2 + ... kn follows Tjf^ — ior q > n. 
Note that due to the coefficients obey 

T,l°) e R , T,!"") = T,!"?) * . (27) 

We stress again that the minimal length of the loops 
in a coefficient T^''-' is n/3, such that we expect the 
ri^'' to decrease exponentially with n. 

The evaluation of the Tn'^^ may numerically be 
implemented efficiently by pre-calculating often used 
combinations such a -ff|.i, H±iHq etc. Note that all 
these terms only live on the single time slice A'^^^ 
Thus the storage requirements are modest and the 
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matrix multiplications and traces needed for the Tn 
run only over the indices in a single time slice. 

In Fig. Owe show the absolute values \Tn'^ \ for var- 
ious q ~ 0, 1...7 as functions of n. We compare two 
values of the mass parameter, using m = 0.05 in the 
top plot and m — 0.20 at the bottom. To illustrate 
the fluctuation of the absolute values we superim- 
pose the data points from three different gauge con- 
figurations. The plots nicely display the expected 
exponential decay with increasing n, with the data 
points for different q and different configurations all 
being close to a straight line. It is interesting to note, 
that although between the top and the bottom plots 
the mass was increased by a factor of 4, the slope 
of this line changes only very little. This is not a 
surprise, since the mass m enters through powers of 
K ~ 1/ (4 -I- m) which for the range of masses consid- 
ered in the plots is only a slowly varying function. 

Also for the phase of the Tn^ we find only a very 
small mass dependence. The weak mass dependence 
of the Tn^ then in turn implies the weak mass de- 
pendence of W which we have observed in Fig. [2] and 
Table I. 



VI. APPLICATIONS OF THE 
FACTORIZATION FORMULA 

The winding decomposition formula can be ap- 
plied to the determination of the canonical parti- 
tion functions for a fixed quark or baryon number. 
The grand canonical partition function Zqc{^j) at 
chemical potential is obtained by integrating the 
determinant over the gauge fields, 

Zacipi) = J P[;7]e-^«[^1 det[D{fi)] , (28) 

where Sg[U] is the gauge action, and for simplicity 
we discuss the case of only one quark fiavor. 

The grand canonical partition function may be 
written as a sum over canonical partition functions 
Z'-'^) with a fixed quark number Q, 

Q B 

where in the second step we have made explicit the 
fact that due to the center symmetry only the canon- 
ical partition functions Z^^^^ , where the quark num- 
ber is a multiple of 3, arc non- vanishing. B is the 
baryon number. We remark, that for expectation 
values also quark sectors where Q is not a multiple 
of 3 contribute, and furthermore stress that the cen- 
ter symmetry may be broken at the phase transition. 

Our winding decomposition formula may now be 
used to determine the canonical partition functions 
in various effective ways. 



A. Projected determinants from the Fourier 
transformation method 

The canonical partition functions may, e.g., be de- 
termined from a Fourier transformation with respect 
to an imaginary chemical potential /i = j<^//3: 

4^^^ - j P[;7]e-^-[^li^('3) , (30) 
where we define the projected determinants Z?^'^\ 

jj(Q) ^ f ^e-"^'^ dei[D(p = t^/p)] . (31) 
J-TT 27r 

Although the determination of the Z)''^) through 
Fourier transformation assumes the elegant form of 
Eq. (j3ip . it is expensive in an actual implementa- 
tion. The (f integral has to be done numerically and 
experience shows that a large number of tp-values 
is needed for an accurate approximation of the (p- 
integral. Thus the determinant det[-D(^ = i(p/P)] 
needs to be evaluated many times for different (p. 

With the factorization formula this procedure may 
be accelerated considerably. Using (jlOp , (fT5|) we find 

DiQ) = ^f\,r'Q^ (32) 

X det[l - i7o - e"^H+i - e-"^H^i]. 

Also here a determinant is integrated over Lp and 
thus has to be evaluated at many values of ^p. How- 
ever, the determinant in ([32]) is for a considerably 
smaller matrix, since Hq,H±i live on only a sin- 
gle timeslice. The strategy thus is to pre-calculate 
Ho,H±i and Aq and to work with p2|) instead of 
([31]) . The dimension of the problem is reduced by a 
factor of (3, the number of timeslices. Since the cost 
for the exact evaluation of a determinant is propor- 
tional to the third power of the matrix size, one gains 
a factor of P^, which for typical lattice sizes is an im- 
provement of 0(100). The increased efficiency due 
to the dimensional reduction allows for a consider- 
ably larger number of values ip used in the numerical 
integration, such that the accuracy of the results for 
the projected determinants improves. 

B. Projected determinants from a power series 
expansion 

In addition to the Fourier transformation tech- 
nique, the equations derived in Section III allow also 
for an efficient determination of the projected deter- 
minants via a power series expansion. From 
(|28)) - ([30| follows the fugacity expansion 

det[i^(^)] = e'^QI^D'^'^^ , (33) 
Q 
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which shows that the D{Q) may be obtained as the 
coefficients of a series expansion of det[D{^)]. 

Using the notation of the last two sections we can 
write the determinant as 

det[D{fi)] = AoW^°^ (34) 

OO 

X Y[ exp (-e^«'^r(«) - e~^«^T(«)*) . 

9=1 

Expanding the exponential functions and ordering 
the terms with respect to the powers of e^^^ we can 
read off the projected determinants D'-'^^ 

(CO OO \ 

9=1 9=1 / 



D^^^ = AaWo -T 



(2) 



^ j.(9+2)y(g)* 
9=1 



et cetera 



(35) 

Again we can obtain the for negative Q via 

complex conjugation, D^"'^^ = D^^^*. Note that a 
term T^"?^ consists of loops that wind at least n — q 
times and thus according to the numerical analysis 
higher terms in (|35p are exponentially suppressed, 
and the sums may be truncated after a few lead- 
ing terms. We remark that the expansion which we 
have denoted here for a single flavor may easily be 
generated to an arbitrary number of flavors. 




FIG. 4: Absolute value (top plot) and phase (j> (bottom) 
of the determinant det[_D] as a function of /i. For an in- 
dividual gauge configuration and three different masses 
we compare the exact result (full curves) to the approx- 
imation discussed in Section VI. C. for A = 1 (dashed 
curves) . 



C. A simple truncation scheme 

Having worked out the evaluation of the contribu- 
tions in the individual winding sectors we now ad- 
dress the question of a suitable truncation scheme. 
Obviously the r(«) are given by infinite sums (see 
Eq. (1221)) which cannot be summed in closed form. 
However, the exponential decrease of the individual 
terms Tn'^\ which was illustrated in Fig.[3l suggests 
that a truncation of the sums might lead to a reason- 
able approximation scheme. We experimented with 
truncated contributions T^^^ defined as 

l+A 

fio) ^ ^ Tio) (36) 

n=l 

kl+A 

f = T^') for q^O. 

n=\q\ 

The parameter A, which is a positive integer, con- 
trols the approximation and in the limit A — > oo the 
exact result is obtained. Since the T^'^^ decay ex- 
ponentially in n, it may be expected that already a 



small value of A might give a good approximation. 
In fact Fig. [3] shows that for this case subsequent 
terms are suppressed by a factor of O(IOO). 

A particularly simple expansion is obtained for 
A = 1. In this case one finds 

f (0) = Tr [Ho] + i Tr [{Hof] + Tr [H+iH^{\ , (37) 

f = i Tr [{H+,r] + Tr [(i/+i)^i/o] , {q > 0). 

Each term (except for q = 0) consists of only two 
traces. Furthermore the arguments of the traces 
have a particularly simple form which allows for a 
efficient recursive implementation of the needed ma- 
trix products. 

In order to test the quality of the truncation at 
A = 1, we evaluated the determinant according to 
(pi]) , replaced the T^"?) by the f («) as given in ([ST]) , 
and took into account all terms with g = 0, 1 ... 6. 

In Fig. m we compare the outcome of the truncated 
determinant (dashed curves) to the exact result (full 
curves). The comparison is done on a single gauge 
configuration for three different masses. The plot 
shows that already the A = 1 approximation works 
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surprisingly well. For most of the /x-range the ap- 
proximation falls exactly on top of the exact curve 
and only when reaches the cutoff at fi = I, devi- 
ations of the approximation (dashed curve) become 
visible. The onset of these deviations may be pushed 
to even higher /i by including larger q in ()34|) . while 
still sticking to the A = 1 approximation. For the 
quenched example studied here this is not necessary, 
since the deviation occurs only close to the cutoff. 

We repeated this analysis for several gauge config- 
urations, and found the same quality of approxima- 
tion throughout. Thus we conclude that the trunca- 
tion p6|) together with the series expansion of Sec- 
tion VLB. provides a powerful method for comput- 
ing the projected determinants Z?^'^) . 

VII. CONCLUDING REMARKS 

In this paper we have revisited the problem of 
finding the fermion determinant in a sector with a 
fixed quark number, i.e., projecting to a fixed wind- 
ing number of the closed loops the determinant con- 
sists of. In our approach we use a partition of the 
lattice into four domains and the determinant is ex- 
pressed in terms of sub-determinants and propaga- 
tors on the domains. The determinant can be orga- 
nized in such a way that the winding numbers of the 
individual contributions can be accessed efficiently. 

Our result allows for two different new approaches 
for computing the projected determinant at a fixed 
winding number: The factorization of the determi- 
nant dimensionally reduces the piece which depends 



on the chemical potential to a single time slice. Thus 
the evaluation of the determinant for different values 
of the imaginary potential f becomes considerably 
cheaper, and many more (^-points can be used in 
the numerical Fourier transformation to the deter- 
minant at fixed quark number. A second more direct 
approach is a power series expansion of the deter- 
minant, where the terms at fixed winding number 
appear as expansion coefficients. We propose a suit- 
able truncation scheme for this expansion which we 
analyzed numerically. 

The method of a partition into domains may 
be used also for other applications. In particular 
it exists also for propagators and allows for the 
analysis of the winding numbers of the paths the 
propagators consist of. This might be useful for the 
recently proposed dual chiral condensate [1], where 
the Fourier transform approach is used to define 
observables that are sensitive to chiral symmetry 
breaking as well as confinement. The domain 
factorization of therm propagator might also be 
combined with a pseudofermion representation 
of the fermion determinant and so open another 
approach to lattice QCD with chemical potential. 
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